Flow reversal in a simple dynamical model of turbulence 
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In this paper, we study a simple hydrodynamical model showing abrupt flow reversals at random 
times. For a suitable range of parameters, we show that the dynamics of flow reversal is accurately 
described by stochastic differential equations, where the noise represents the effect of turbulence. 
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It has been recently reported Q that abrupt flow re- 
versal takes place at large Rayleigh number in thermal 
convection. Beside thermal convection, flow reversal has 
been also observed in laboratory experiment of two di- 
mensional turbulence |2| and in the magnetic polarity of 
the earth More generally, there are many "turbu- 
lent" flows for which transitions between different states 
have been investigated, namely within the theory of mul- 
tiple equilibria for atmospheric flows £| and of long time 
climatic changes 

In most cases the major question to be answered by 
experiments or observations concerns the mechanism re- 
sponsible for the transition. This question is highly 
non trivial whenever the average persistent time < r > 
around different states is much longer than any charac- 
teristic times describing the dynamical behavior of the 
system. 

There are two main interpretations which have been 
suggested so far. The first one assumes that turbulence 
is more or less a "noise" applied to the order parameter 
-0 which describes the system (i.e. the wind in the case of 
RB convection or the temperature in the case of climatic 
change) . Then, by defining ip such that the two observed 
states are iipo, the equation for ip is given by: 

tb 2 

dtp = [aip(l - -^)]dt + s/^dWit) (1) 

n 

where W is a Wiener process S correlated in time while 
a" 1 is the characteristic time scale of the instability at 
tp = 0. One can show that transitions between the two 
stable states iipo occur at random times r with an aver- 
age time < r > given by 



< r >r 



s/2a 



exp(aip 2 /(2a)) 



(2) 



Hereafter, following the language of stochastic differential 
equations, the random time r will be referred to as "exit 
time" . Note that equation (fjj implies, for small cr, that 
the average exit time < r > is much longer than the 
deterministic " fluidodynamical" time a~b Moreover, by 
employing the theory of stochastic differential equations 
0, one can compute the probability distribution of the 
exit time r, which, for small cr, is given by: 



P(t) =< t > 1 exp(-r/ < r >) 



(3) 



repeated small noise perturbations of the same "sign" 
which are acting against the deterministic "force", i.e. 
there is no specific mechanism introduced by the small 
scale turbulence (parametrized by the noise) and transi- 
tions can be explained in terms of large deviation theory. 

One of the major criticism against the above interpre- 
tation is that the noise is by itself the effect of turbulence 
and, in most cases, there is no time scale separation be- 
tween the dynamic fluctuations of ip and turbulent fluc- 
tuations. Thus one cannot assume that "the noise" is 
a "fast" perturbutaion with respect to the dynamics of 
ip and, as a consequence, equation JIJ cannot be justi- 
fied. As an alternative way, one should look for a specific 
"fluidodynamical" large scale mechanism which can ex- 
plain the observed transitions. For instance, for the wind 
reversal in thermal convection, there has been a recent 
proposal which explains transitions as the result of 
plume dynamics. 

In this letter we want to understand whether and how 
equation (Q can be justified, at least in the simplest pos- 
sible model of a "turbulent" flow. For this purpose we 
shall consider an 'energy cascade' model i.e. a shell model 
aimed at reproducing few of the relevant characteristic 
features of the statistical properties of the Navier Stokes 
equations In a shell models, the basic variables 

describing the 'velocity field' at scale r n = 2~ n r = k^ 1 , 
is a complex number u n satisfying a suitable set of non 
linear equations. There are many version of shell models 
which have been introduced in literature. Here we choose 
the one referred to as Sabra shell model 



= ik n [aku* n+1 u n+2 

dt 



-2U„_1 



+ fn 



(4) 



where A = 2, a = 1 and c = —(1 + b) and /„ is an exter- 
nal forcing. Let us remark that the statistical properties 
of intermittent fluctuations, computed either using u n or 
the instantaneous rate of energy dissipation, are in close 
qualitative and quantitative agreement with those mea- 
sured in laboratory experiments, for homogeneous and 
isotropic turbulence |10|. 

The basic idea of our approach is to assume that 
U r = real{u\) describes the one dimensional unstable 
manifold arising by a (large scale) pitchfork bifurcation. 
Consequently we change the equation for u\ as follows: 



If the above scenario is believed to be correct, then 
the transitions between the two states iipo are due to 



— = $ + - ukfui 

dt Un 



(5) 
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= ik n [aAu* l+1 u n+2 + bu^_ x u n+ i - cA 1 u n -2U n - 1 ] 
-vk 2 n u n + /„ . (n > 1) (6) 

where $ = ifciaAu2M3- Let us comment equation (JSJ). In 
most cases, a pitchfork bifurcation, as described by equa- 
tion is observed in real fluidodyanamical flows with 
respect to the external forcing or the Reynolds number. 
In our case we are assuming that the unstable manifold 
is coupled to smaller scales by the term ikiaAu^u^. For 
small v, the two states U\ = ±Mq becomes unstable and 
a turbulent regime is observed. In the following we will 
think of equation JSJ as a realistic, although approxi- 
mate, equation describing a simplified turbulent "flow" 
superimposed to a large scale instability. As one can see, 
no external noise is introduced in the system. 

We remark that equations 1516(1 are symmetric un- 
der the transformation Ui — > —U\. More precisely, one 
can show that by changing (u 3m+ i, u 3m+2 , u 3m+3 ) with 
(-u 3m+ i,-u 3m+2 , u 3m+3 ) ( m = 0, 1, 2, ..), the equations 
of motion are invariant. Using dimensionless variables 
W = u n /uo, K n = k n L and t' — fit (L = fcf 1 ), one gets: 

ItF = ^ K l aAW n+l W n+2+W 1 (l-W?)-j^KlW 1 

Equation Q tells us that the dynamical behavior of 
u n depends on two dimensionless number, namely the 
Reynolds number Re = uoL/v and the number B = 
Uq/(u.L). We will investigate the statistical properties of 
eq. (l- r )l(H for Re — > oo and for different values of B. In 
particular we fix fi — 1 and uo real, while the parameters 
of the model are a = 1, b = —0.4, c = 0.6. 
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FIG. 2: Same as in figure Q for ito = 2. 

abrupt reversals of U r are observed at apparently random 
times in both cases. The most important feature of figure 
Q and J2J) is that the characteristic correlation time of 
U r is of order 1 (i.e. it is of order L/u ~ 1), much 
smaller than the exit time r. The behavior shown in 
figures ill 121) does not change by increasing the Reynolds 
number. A more refined numerical simulation shows that 
the "random" exit times r are distribute according to 
eq.© with < t >~ 600 and < r 65 for B = 1 and 
B = 2 respectively. In figure we show logP(r) versus 
r for the case B = 2, where the line shown in the figure 
represent the quantity exp{— r/65). 
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FIG. 1: The velocity U r plotted as a function of time as ob- 
tained by numerical simulation of equation© for uo = 1, 
H = 1 and v = 10" 6 . 

In figure we show U r as a function of time for v = 
10~ 6 and B — 1 (uq — 1), while in figure we show 
the same variable for B = 2 (uq = 2). As one can see 



FIG. 3: log[P(r)\ as a function of r as obtained by numer- 
ical simulations of @ for uo = 1. According to the theory 
of stochastic differential equations, P(r) is exp(—r/ < r >) 
where < r > is the average exit times. The line in the fig- 
ure shows exp(— r/65) which is an extremely good fit to the 
observed -P(t). 

For larger values of mq, the average exit time < r > 
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becomes smaller and eventually, for uq ~ 10 it becomes of 
order 1 (see also the discussion below). Our interest will 
focus for values of u$ in the range [1, 2] where < r > is at 
least two order of magnitude longer than and L/uq, 
the two relevant deterministic time scales of equation (J5J). 

A closer look at figures JU and (J2J reveals that the two 
states of U r are not ±uq. For B = 1 the maxima ±Um 
in the probability distribution of U r are located at Um = 
0.84, while for B = 2 we find Um = 1-5, represented as 
horizontal lines in figure J2J. 

In order to explain these results, let us consider more 
carefully the physical meaning of $ in equation (jSJ) . The 
quantity real($>u\) is the amount of energy transferred 
by mode ui to smaller scales, i.e. U2 and U3. On the 
average, we know that real < (^u*) >= — e < 0, where e 
is the average rate of energy dissipation. Thus, as a first 
approximation, we can assume that 



$ = -(3ui 



(8) 



where < <j>'u\ >= and (3 > 0. By multiplying both side 
by u*, taking the time average, we can compute (3 as: 



/3 = real(< $ul >)/ < |iii| 2 > 



(9) 



Thus we must expect that two maxima in the probability 
distribution of U r should corresponds to the solution of 
the equation: 



(/*■ 



< M > 



)u r -^ul = 0. 



(10) 



Equation IjlUI) tells us two important informations. First, 
the "states" ±£7m> (between which abrupt transitions are 
observed) are not stationary solutions of the determin- 
istic equations JSJ, rather the "states" should be con- 
sidered as statistically stationary states of the system. 
Second, we must expect Um < «o a s far as an energy 
cascade, e > 0, is produced. We have carefully checked 
the validity of equation ifTUI) for u in the range [1, 2]. We 
have computed (3 from numerical simulations according 
to equation |(?JJ|. It turns out that the numerical val- 
ues of (3 are extremely well represented by the function 
P(uo) = 0.15+0.14uq. Next we have solved equation fTT)|) 
for each value of uo and (3(uq) getting the line represented 
in figure J3J . Finally we have computed the value of Um 
using the the probability distribution of U r as obtained 
by the numerical simulations of the model. The results 
are shown as symbols in figure Q: the agreement is ex- 
cellent. 

At the light of the above results, it is tempting to argue 
that the behavior of © is consistent with the stochastic 
differential equation: 

u 2 

dm = [U - 0(u o ))ui + fj,ui-^]dt + y/adW(t) (11) 

for a suitable value of the noise variance a. Using @ and 
|J5J), the quantity yfadW(t) represents the deterministic 
term (f>'(t)dt. In order to validate (flip, we need to check 




FIG. 4: Computing the statistically stationary solutions for 
U r - The line corresponds to uo y/fx — /3(wo), where P(uo) is 
defined by equation The symbols represent the value of 
the maxima Um of the probability density function of U r ob- 
tained by the numerical simulations of JSJ for different values 
of Up . 



whether the observed fluctuations of u\ are in agreement 
with the average exit time < t >. More specifically, 
using l|ll|) and we have: 



< t >= 



V2( M - 0) 



exp({fi - I3) 2 ul/2^a) 



(12) 



By the numerical value of < t > and (3, we can com- 
pute the amplitude of the noise tr, hereafter denoted by 
<7 T , needed to explain the observed average exit times. 
The crucial point is whether the observed fluctuations 
of U r near one the two states are "compatible" with the 
noise amplitude a T . To answer this question, let 5U be 
small deviation around the statistically stationary states. 
Using lfTT|) we can compute < 8U) 2 > as: 



< (SU) 2 >-- 



4(/i - (3) 



(13) 



By using the numerical simulations, we can estimated 
< (SU) 2 > around the statistically stationary states. Fi- 
nally, using (|13|l . we can estimated the "noise" variance 
a, hereafter referred to as asu, which explains the ob- 
served fluctuations of SU. For equation (|11|) to represent 
a good approximation of the full non linear deterministic 
system, we must obtained 



0~5U 



(14) 



In figure JSJ, we plot a T (line) and crsu (crosses) for 
different values of uo in the range [1,2]. As one can see, 
(|14fl is verified with very good accuracy, with at most 10 
per cent difference for uq = 2. We remark that the result 
shown in figure |J5} represent a rather severe test on the 
validity of equation 111|) . 
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FIG. 5: Plot of a&u (symbols) and a T (line) for different values 
of uq.(t t is the noise computed by equation 1121 while asu 
is the noise computed from the fluctuations of U r near the 
statistically stationary states, see equation l|13|l . 

We want finally to comment on the behavior of the 
model in the limit of large uq or, equivalently, in the 
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limit of large B. For large value of uq we reach the con- 
dition /3(u ) ~ [i. Actually, arguments based on energy 
balance and numerical simulations clearly show that the 
maximum value of is 1. In this case, equation (|1(J|) gives 
Um ~ 0, i.e. the two "statistically stationary states" dis- 
appear. It follows that, for large B, one cannot speak of 
"abrupt flow reversal" . 

All the results, presented so far, strongly indicate that 
equation is a very good candidate to explain the 
observed abrupt "flow reversal" shown in figures and 
(J2J), i.e. "flow reversal" is explained within the framework 
of large deviation theory of stochastic differential equa- 
tions. Therefore, we argue that some of the experimental 
results mentioned in the introduction, can be investigated 
and explained in the framework of stochastic differential 
equations even if no time scale separation exists between 
the "wind" fluctuations and turbulence. 
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